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of. .a. Sj;h.c?ic.eA  Explosion  in  Sea  Watpg* 


ty 

D.  K.  Y.  Ai^and  H.  Holt2 


Abstract. 


This  report,  which  Is  a  sequel  to  papers  by  Holt  (1955 
and  1956),  describes  the  first  stages  of  the  calculation 
of  the  growth  of  a  spherical  explosion  in  sea  water,  due 
to  a  charge  of  PETN  initiated  at  its  center ,  Working  in 
the  time  distance  plane,  the  analysis  of  Holt  (1955  and 
1956)  is  used  to  construct  a  "plus"  characteristic  travers¬ 
ing  the  whole  field  of  disturbance  near  the  origin  of 
blast  and  the  "minus"  characteristic  which  is  the  outer 
boundary  of  the  detonation  regions  These  initial  data 
are  then  used  to  start  a  calculation  of  the  further  growth 
of  the  blast  field  by  the  numerical  method  of  characteris¬ 
tics,  This  is  programmed  on  an  IBM  CPC  computer,  A 
procedure  is  established  which  is  satisfactory  in  principle 
but  unsuitable  for  repeated  application  on  this  particular 
computer  owing  to  inadequate  storage  facilities.  As  a 
consequence  only  the  first  stages  of  the  calculation  are 
carried  out  now  and 'the  work  will  be  completed  later  on 
an  IBM  7cU  computer. 


*  The  results  presented  in  this  paper  were  obtained  in  the 
course  of  research  sponsored  by  the  Office  of  Naval  Research 
under  Contract  Nonr-562(07) 

1,  Research  Assistant  in  Applied  Mathematics 

2,  Assistant  Professor  of  Applied  Mathematics 
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I,  Introduction. 

This  report  is  a  continuation  of  work  by  Holt  (1955  and 
1956)  or  the  initial  behavior  of  a  spherical  explosion.  A 
beginning  is  made  with  the  numerical  calculation  of  the  growth 
of  a  spherical  explosion  in  sea  water,  due  to  a  charge  oi  PETN 
initiated  dt  its  center.  The  basic  data  are  taken  from  Holt 
(1956),  where  series  expansions  are  developed  to  determine  the 
field  of  disturbance  near  the  origin  of  blast  (the  point  0* 
in  the  time  distance  piano  shown  in  Figure  1) . 

Three  stages  preliminary  to  this  calculation  are  completed 
in  the  present  reports  Firstly,  the  series  expansions  are  used 
to  calculate  initial  values  of  physical  variables  on  a  "plus” 
characteristic  traversing  the  whole  disturbance  field  near  O'. 
Secondly,  the  ordinary  differential  equations  satisfied  in  the 
detonation  region  are  integrated,  using  Jones'  equation  of 
state  for  PETN}  the  integration  leads  to  the  contst ruction  "df 
the  minus  characteristic  dividing  the  detonation  region  from 
the  gas  expansion  region.  Thirdly,  the  data  on  these  two 
characteristic  lines  are  used  to  initiate  a  calculation  of  the 
remaining  field  of  disturbance  by  the  numerical  method  of 
characteristics.  This  work  is  programmed  on  an  IBM  CPC  computer 
attached  to  the  Division.  The  program  is  satisfactory  in 
principle  but  requires  a  larger  storage  capacity  than  that 
available  on  this  machine.  Accordingly  only  the  first  stage 
of  the  numerical  calculation  has  been  carried  out  so  far  and 
the  program  will  be  completed  on  an  IBM  ?CU  computer  which  has 
the  required  storage  facilities. 
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The  physical  phenomena  and  t.he  theoretical  analysis  of 
the  problem  are  given  below  briefly#  The  description  is 
derived  from  the  fuller  account  given  in  Berry  and  Holt  (195*+)  • 

When  a  spherical  charge  of  orthodox  explosive  is  initiated 
at  its  center  it  is  observed  that  the  main  blast  wave  is 
followed  by  a  second  blast  wave  which  is  very  weak  at  first 
but  attains  an  appreciable  strength  later#  The  cause  of  this 
second  wave  is  the  slightly  excessive  expansion  of  the  explosive 
gases  behind  the  detonation  front  which,  if  unchecked,  would 
make  the  gas  pressure  at  its  outer  boundary  too  low  in  comparison 
with  the  adjacent  fluid  pressure  behind  the  powerful  main  blast 
wave#  The  second  blast  wave,  which  always  starts  Inside  the 
explosive  gas  regionr  has  the  effect  of  increasing  the  gas 
pressure  sufficiently  to  ensure  continuity  at  the  fluid-gas 
boundary# 

Theoretical  analyses  of  the  existence  of  the  second  blast 
wave  have  been  done  by  Whifcham  (1950),  Wee ken  (1951),  Berry 
and  Holt  (1955:),  Berry,  Butler  and  Holt  (195*0  and  Holt  (1955) 
based  on  the  investigations  made  independently  by  Taylor  (1950) 
and  Doring  (Dcring  and  Burkhardt  (19^6) )«  The  early  develop¬ 
ment  of  the  disturbance  from  a  typical  spherical  explosion  is 
illustrated  in  Figure  1,  which  shows  trajectories  in  the  t,r 
plane,  where  t  is  the  time  measured  from  the  instant  of  initia¬ 
tion  and  r  is  the  radial  distance  from  the  center  of  the  explo¬ 
sion#  It  begins  with  a  detonation  phase,  during  which  a  strong 
detonation  wave  0fD  travels  outwards  from  the  center  and  reacts 
on  the  solid  explosive  to  produce  highly  compressed  gases# 
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This  is  the  main  blast  wave,  which  is  followed  by  a  region  of 
disturbed  fluid  denoted  by  c„  At  the  same  time  the  gases 
released  by  the  explosion  expand  rapidly  through  a  centered 
expansion  wave  A0!B,  a  region  denoted  by  e«  At  the  head  of 
this  expansion  wave  is  another  gas  region  of  initially  uniform 
flow  BO'C  denoted  by  g,  this  is  adjacent  to  the  compressed 
fluid  region  c* 

In  the  stag9  following  the  detonation  process  character¬ 
istic  curves  play  an  important  part,  We  define  a  "minus 
characteristic"  to  be  a  line  in  the  t,r  plane,  the  slope  of 
which  is  everywhere  equal  to  the  difference  between  the  local 
velocity  of  fluid  and  the  velocity  of  sound.  Correspondingly, 
a  slope  of  a  "plus  characteristic"  is  everywhere  equal  to  the 
sum  of  these  velocities.  Throughout-  this  report  we  are  dealing 
only  with  these  two  types  of  characteristic.  In  water  the 
entropy  change  is  considered  to  be  negligible,  therefore  the 
third  type,  a  streamline,  which  is  a  characteristic  for  the 
entropy,  does  not  arise 0  Entropy  changes  are  significant  in 
region  g,  but  not  at  the  small  distances  from  0*  considered 
here.  In  the  expansion  region  a  fan  cf  minus  characteristics 
radiates  from  0*.  This  whole  region  is  isentrooic  and  can  be 
described  by  tx<;o  dependent  variables  u  and  a.  Region  g  is 
bounded  by  a  minus  characteristic  line  on  one  side  and  a  stream¬ 
line  on  the  other.  This  region  Is  isentropic  only  to  the 
extent  to  which  its  boundary  with  region  e  remains  a  minus 
characteristic.  In  the  compressed  water  region  the  flow  is 
entirely  isentropic.  In  both  regions  three  dependent  '/ariables 
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u,  p  and  a  are  used* 

In  the  neighborhood  of  O’,  equations  of  motion  are  solved 
in  series  expansions  to  describe  the  state  of  motion#  Further 
out  in  the  field  the  method  of  characteristics  is  used#  Weoken 
(1951)  employs  series  expansions  in  the  neighborhood  of  0*  ofrjty 
up  to  terms  or  order  (t  -  t0)*,  where  (t0,  r0)  are  the  coor¬ 
dinates  of  O’,  relying  on  numerical  work  to  continue  the  field 
further  out#  This  is  inadequate  to  gain  a  full  understanding 
of  the  initial  field  of  disturbance,  and,  in  particular,  to 
prove  conclusively  that  a  second  shock  must  always  develop, 
the  series  must  be  taken  to  higher  order  terms #  Berry  and 
Holt  (1954)  made  an  investigation  of  the  mathematical  nature 
of  the  singularity  in  full#  The  flow  is  expanded  in  half¬ 
powers  of  the  distance  from  O'  in  the  t,r  plane  with  coefficients 
depending  on  the  angular  coordinates#  The  expansions  are  not 
valid  in  the  neighborhood  of  the  minus  characteristic  through 
O’  in  region  g.  Near  this  line  the  expansions  are  modified  by 
the  extension  of  a  technique  due  to  Lighthill  (1949)#  The 
method  developed  by  Berry  and  Holt  (1954)  was  generalized  by 
Holt  (1955)  so  that  a  soherical  explosion  in  water  may  be  con¬ 
sidered#  The  result  of  their  work  is  summarized  in  II# 

In  this  report  calculations  are  made  to  a  spherical  charge 
in  water  based  on  the  method  derived  by  Holt  (1955)#  The 
series  expansions  are  carried  out  separately  in  each  of  the 
regions  and  coefficients  are  computed  according  to  the  boundary 
conditions  and  equations  of  state  in  Holt  (1956)#  One  member 
each  of  the  plus  and  minus  characteristic  lines  are  constructed# 
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The  plus  characteristic  starts  at  a  point  on  the  detonation 
front  very  near  O’  and  goes  through  the  whole  disturbed  region,. 
The  minus  characteristic  goes  radially  from  Of  along  the  expan¬ 
sion..-  detonation  boundary,,  With  the  positions  and  conditions 
of  these  two  characteristics  known,  further  out  in  the  region 
other  members  of  the  characteristic  families  may  be  constructed 
in  a  step  by  step  process e  Numerical  integration  and  an  itera¬ 
tive  method  are  employed  and  programmed  on  an  IBM-C PC  computer,, 
It  is  found  that  even  with  the  aid  of  the  IBM  computer  the 
problem  requires  a  considerable  amount  of  time  to  solve  com¬ 
pletely,  However,  all  results  found  are  tabulated  so  that  we 
are  in  a  ready  position  to  extend  the  work* 

11  •  1*3, aasjpn  in  series  jmzJ&e,  origin  of  the  main  blast  wave. 
The  equations  of  momentum  and  continuity  In  unsteady 
spherical  symmetric  flow  are,  respectively 


3p 

at 


+  P 


1  Oja  „ 

p  3r  " 


0 


u  M  =  _  ?  aa 

or  r 


Throughout  the  analysis,  all  variables  are  made  non-dimensional, 
The  radial  distance  r  and  the  time  t  are  divided  by  their 
respective  values  at  the  surface  of  the  sphere,  the  pressure  p 
and  the  density  p  by  their  respective  values  at  the  detonation 
front  and  the  fluid  velocity  u  and  velocity  of  sound  a  by  the 
detonation  velocity. 

New  variables  £  and  8  are  introduced  and  defined  by 
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r  =  1  -  ?2sin  8 
t  =  i  -  ^2cos  e 


£2  and  9  are  polar  coordinates  in  the  t,r  plane  based  on  O’, 

2 

with  9  equal  to  the  angle  measured  below  0rM  and  £  the  radial 
distance  from  0*9 

The  form  of  expansion  to  be  assumed  for  each  dependent 

variable  is  determined  mainly  by  the  known  behavior  in  the 

detonation  region.  Using  the  results  of  Taylor  (1950)  and 

Doring  (Doring  and  Burkhardt  (19*+6))  the  variables  are  analytic 

functions  of  C  where  C  =  9.~  9Q?  9)^  , 

(1  -  £rcos  9)t 

Berry  and  Holt  (195*+)  assume  series  expansions  for  u,  a 
and  p  in  the  remaining  regions,  of  the  form 


u  =  u0(@)  +  ux(eK  +  ^(9)^  +  »»9 

a  =  a0(9)  +  axCeK  +  a2(9K2  +  ... 

p  =  po(0)  +  Px(eH  +  P2^H2  +  ••• 

The  expansions  for  each  region  are  different,  and  are 
distinguished,  one  group  from  another,  by  using  an  appropriate 
suffix. 

The  Isent.ronic  region. 

In  the  isentropic  region  e,  after  introducing  the 
auxiliary  variable  of  state 


where 


d 

k 


=  l 

k 


is 

pa 


P 
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the  equations  of  momentum  and  continuity,  when  expressed  in 
terms  of  the  independent  variables  and  9,  may  be  written  in 
the  form 

cos  9  +  (u  +  a) sin  9?(4£  + 
l  J  3S  9S 

,2 


(2.1) 


-  2jsin  9  -  (u  +  aicosl©  C|{j  ±  |^j)  +  - — — -  =  0 


l 


1  -  S  sin  9 


In  equations  (2.1)  after  substituting  expansions  of  the 
dependent  variables  of  the  form 

ae  =  uoe(9)  +  aie(0)*-  +  u2e(e)^2  + 

and  equating  corresponding  coefficients  of  £ ,  one  group  of 
equations  are  obtained  for  each  order  of  coefficients.  In  this 
paper,  the  n’th  order  coefficient  is  referred  to  as  the  coeffi¬ 
cient  of  £n. 

Terms  independent  of  K  in  equations  (2.1)  give  the  rela¬ 
tion  between  zero-order  coefficients 

uoe  *  aoe  =  tan  9  (2.2) 

uoe  +  doe  =  constant  (2.3) 

The  terms  of  order  K  in  (2.1)  give  the  following  rela¬ 
tions  between  the  first-order  coefficients* 

(d-^g  -  ^2.0^  2^aXe  “  ale^  ^ule  "*  ”  0  (2.h) 

|oo3  a  +  (uoe  +  aoe)sin  9  }(ule  +  ale) 

-  2[sin  9  -  (uoe  +  a„e)ooS  a](u£e  +  8]9>  =  0 


(2o?) 
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where 


u*  =  ,  etc, 

oe  7 


d© 


p 

The  relation  between  the  coefficients  of  £  derived  from  (2.1) 
are 

2(u0  -  d  )  +  2(u*  -  d  *  )  ( u  -  aD  )cos2© 

2e  2e  oe  oe  2e  2e 

+  (ule  ~  ale)(ule  “°le)sin  @  cos  6 


+  2(ule  -  ale)(u'e  -  dje)cos  9 


+  4aoeaoecos  9=0 


(2.6) 


2 -(cos  ©  +  (uQe  +  aQe)sin  © 


j(u 


2e  +  d2e) 


-  2isin  8  -  (a  +  a  )cos  9Hu'  +  d  *  ) 
oe  r.->  J  2e  2e 

+  (ule  +  +  °Le  )sin  9  +  2(uie  +  <’le)o°s  9) 

^Uoeaoe  -  ®  (2.7) 


Both  variables  d  and  a  are  related  to  p  by  the 

oe  oe  oe 

equation  of  state  so  that  pQe  is  defined  implicitly  by  the 
equation 

°oe(Poe)  +  aoe(Poe)  =  const>  “  tan  ®  <2*8) 

The  equation  of  state  in  region  e  may  be  written 

P  =  p(p) 


which  becomes,  on  expansion 
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P  +  p,  £  +  p,  [ 
oe  les 


le 


=  p(p  )  +  ka^  p_  C 
roe  oer  le 


+  (&koep2e  +  *Leple^  ao° 


where 


l°e  "  vdp  oe 
We  then  have  the  relations 


kaoe  =  »  Lfi  =  ^^f^oe 


& 

dp' 


p_  =  ka  pn 
le  oe^le 

p2e  =  kaoep2e  +  ^Leple 


Similarly,  Iby  expanding  the  equation 


&  =  JL- 

dp  kpa 

we  find  that  pl9  =  ka09Poe3l9 


p2e  =  k(aoepoe°2e  +  *°ieCaoePle  +  poeale)] 


Prom  the  equation  ka2  =  we  obtain 

dp 


ale  ~  2kal0  pl© 


a 


2s  2ka 


1—j 


oe 


Lep2e  +  *  ka?. 


le 


Combining  the  previous  results,  we  find  that 

ale  -  Vie 


X  = 


_  If, 


2ka‘ 


where 
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The  above  relation'  is  useful  in  solving  ule,  a-^  and  d 
in  terms  of  zero-order  coefficients,  A  similar  procedure  can 
then  be  applied  to  equations  (2,6)  and  (2,7)  to  find  u2e»  a2e 
d2e*  ac^ual  calculation  of  these  coefficients  will  be 

presented  in  detail  in  III, 


Region  c  is,  in  general,  entirely  non-isentropic  and  region 
g  is  partially  so.  Much  of  the  analysis  of  region  g  can  be 
deduced  as  a  simplification  of  that  in  region  c. 

After  changing  the  independent  variables  r,t  to  £  and  9, 
the  equations  of  momentum  and  continuity  are  in  the  following 
form* 


^{cos  0  +  (u  +  a)sin  (kpa  Ql  +  -|ll) 

(  J 


-2(, 


2|sin  9  -  (u  +  a) cos  9 


tkpa  Is  1  ^ 


*ttaa2p,V 

1  -52sin  9 


=  0 


(2,8) 


In  the  above  equations  we  substitute  expansions  for  each 
dependent  variable  of  the  type 

u  =  uoc  +  ulc(9)5  +  u2c(©U2  +  ..o 

p 

where  u„  Is  a  constant.  Equating  coefficients  of  K  and  K  % 
we  then  have  a  set  of  ordinary  first-order  differential  equa¬ 
tions  for  the  coefficients  of  first;  and  second  order  in  turn. 
Solving  these  equations  we  obtain 
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uio  -  irhr-  Wj*t +  wh 

^  oc  oc 
pic  ~  °i't|i - 

Pie  ■  ^5  <C1+1  '  °i*£> 

ka 


a 


1c 


=  L|«  - .  (e^f  -  Co^t'! 
2k2a3  1  11  2  2  I 

00 


J 


where 


=  sin  8  -  (uoc  +  aoc)cos  9 
^2  =  (uoc  “  aoc)  "  sin  9 


-  (M) 


8p 


2  oc 


and  the  are  integration  constants* 

To  find  coefficients  of  <2,  auxiliary  variables  q  and  r 
are  introduced  and  defined  by 

q  =  kp0aQu  -  p 
r  =  kp0aQu  +  p 

We  obtain  from  the  corresponding  coefficients  6£  5  2  irt< (2e8) 

r?c  *  |2k2u0CaQCP oc  -  <1  +  MCC?  +  c|)j 

^OO  oc  *• 
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u2c  and  p2c  may  "be  solved  from 


q2c  =  kpoaott2c  '  P2= 


"2c  =  kPoaou2e  -  p2o 


The  remaining  second-order  coefficients  pgc  and  agc 
can  be  deduced  from  the  expansion  of  the  equation  of  state 
p  =  p(p,S)  and  ka2  =  3p/0p  at  the  initial  condition#  We  obtain 
the  relations 

tag  =  +  £PP®  -  ka®) 


where 


p  =  <«»)„ 

3p’3  0 


'EM  parUallyrigentpQplc ,  j_e.gj„gfl 

In  region  g,  the  situation  is  quite  similar  to  that  of 
region  c,  The  zero-order  coefficients  are  constants.  The  first- 
order  coefficients  are 


Ulg  =  ^"V"  <0l,’l  '  °2,>2> 

^og  og 

Pig  -  “if  i  -  0s?I 

plg  =  ^2“  (0l'Pl  -  °2<pf> 
og 
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The  second-order  coefficients  are  obtained  through  the 
auxiliary  variables  r  and  q. 


=  kp  a  u_  + 

2g  *  og  og  2g  ^2g 


cos  9 

kp  a 
og  og 


|2k2u 


-  (1  +  \„)(G2  +  G2) 
og  og  og  g  1  2  . 


3  P2 

1-2.3/2 


G,  G0Cpn  r 

3kpoga5gl  g  9i 


-  3d  +  Xg)(^/ 


■2Vi 


q0  -  kp  a  u_,  -p_ 

2g  Hog  og  2g  ^2g 


=  -Sfifi-SL  X(1  +  \J(G2  +  G2)  -  2k2u  a3  92  1 

kp  a  1  g  1  2  og  ogrog  f 

\og  I  J 


Og  og 


where 


be 


<p1  =  sin  9  -  (u  +  a  )cos  9 

X  og  og 

<p0  =  (u  -  a  )cos  9  -sin  9 
2  og  og 

and  the  G*s  are  integration  constants.  Pgg  and  a^g  can 

solved  in  the  same  manner  as  p.  and  a.  a 

2g  2c 

The  expansion  of  r  is  valid  throughout  region  g;  but  the 
expansion  of  q  is  not  valid  along  the  minus  characteristic 
through  0*.  In  order  to  overcome  this  difficulty,  an 
independent  variable  z  is  introduced,  where 


9=6  +  z  +  6lg(z)5  +  90n(.zX  +  *•* 


2g 


6  is  the  initial  value  of  9  on  the  minus  characteristic  line 
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through  O’. 


The  new  expression  of  q2g  in  terms  of  z  may  he  written 


,(z)  = 


+  V(°i +  Gi>  -  a2,iog^gpog) 

j(i  .  K  )(2i)3/2  .  3(1  +  x  ,(21)*] 
3kp  «?_  I  ®  'PS  2  ?2  J 


3kp  ac 
og  og 


3  kp  a2 
*og  og 


tt  cos  6(1  -  Xg)  (2aog)  2sin  2z  cos  z  -  2G^92 


where  9^  and  9^  are  evaluated  for  the  argument  6  +  z,  At  z  a  0 


or  9  =  6 


q  =  ...QQS& — id  +  X  )(G2  +  G  )2  -  2k2u  a^  p2  ~  2G,-90 
2g  kpogaog  [  S  1  2  og  ogKogJ  5^2 


coext lc: 


The  equation  of  state  in  dimensionless  form  in  region  e 


j  =  1,6309  -  l,242lp  +  0„5l65p  +  0.09>+56p“5/6  (3.1) 

From  (3,1)  we  deduce  the  velocity  of  sound 


a  =  (1  4»)* 
k  dp 


=  (1)*  X&m  _  Ox24aiiL± 

k  (1,2421  -  1.033P  +  0,078875p“11'/6)t 


,-5/6 


(3,2) 


and  the  auxiliary  variable  of  state 
P  P 

if  f  ' 

=  k  ^  =  (1,2421  -  l,033p  + 

o  0 


033p  +  0,078875p"11/6)^dp 
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where 


k  =  D2pd/pd  =  5.^09 


Knowing  a  and  o  as  functions  of  p  we  can  next  find  9  as  a 
function  of  the  initial  value  of  p  in  the  expansion  zone,  pQe, 
from  the  relation 


‘WPoq5  +  ‘WPoe5  =  oonstant  "  tan  0  (3*3) 


tfhe  constant  is  immediately  found  from  data  at  the 
detonation  front  at  which  pQ0  =  0.6132,  8  =  2.8876,  °oe(Poe) 

=  1.49186,  aQe  =  0.5787?  constant  =  1.81101. 

All  the  initial  values  of  dependent  variables  in  the 
expansion  zone  are  then  defined  implicitly  as  functions  of  9. 

To  obtain  the  first-order  coefficients,  consider  first 
equation  (2.5),  the  solution  tf  which  gives  us  the  ratio  of 
ule(9)  +  d^Q(0)  to  ule(cO  +  dle(°0.  we  may  write 


Ule(e>+01.(8 


)-/u 


ft 


008  8  +  (ttoe*aoe 


)sin  9] 


le(a)*ie(a)\e 


(sin  6  -  (uoe+aoe)cos  6 

(3’ 


d0 

4) 


a  is  the  initial  value  of  the  boundary  between  e  and  d. 

In  order  to  determine  ul8(o)  +  dl0(a),  we  take  the  result 
of  Dorlng  that  in  the  detonation  region  u  and  d  are  functions 
of  C  and  may  be  written 


2 

U  ~  Uq  U^C  U2C  t  « 0 . 
d  ~  d8  dg^*  .8. 


where  the  coefficients  are  constant  and 
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r2  -  S.  (PlR.fi- =...fiQP  ei 
2 

1-5  cos  9 

At  the  boundary  between  regions  e  and  g 

©  =  a  +  +  $2Z2  +  ... 

Novr,  expand  C  in  Taylor* s  series  about  the  point 


©  =  a 


we  obtain 


,  -  8,  (cos  a  +  sin  a) 

C  =  5(sin  a  -  cos  a) a  +  52  _ _ _ _ _ 

a)* 


+  •< 


2(sin  a  -  cos 


Substitute  into  u  +  d  to  get 


u  +  cj  =  (uQ  +  <Jq)  +  (u.  <J1)(sin  a  -  cos  a)?£ 


{3-,(cos  a  +  sin  a) 

+  [(u2  +  d2)(sin  a  -  cos  a)  +  (u.  +  d  ) - ,]5‘d+.,, 

•*-  r\  r  ~  t _  _  No* 


2 (sin  a  -  cos  a/3 


But  in  region  e  we  have 


ue  =  uoe<9)  +  ule(9)C  +  u2e(9)i-  +  ”• 
°e  '  '5oe(e)  +  cle(8H  +  +  ••• 


Expand  ue  +  d^  about  ©  -  a  and  compare  coefficients  of  same 
powers  of  5  with  those  in  u  +  d „  We  obtain 

uo  +  do  =  +  <*—(«) 


oe 


oe 


(Ui  +  ) (sin  a  -  cos  a)2  =  n^a)  +  d^^a) 

Considering  equation  (2,4)  and  taking  into  account  the 
previous  results  we  obtain 
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uis(9)  +  %(e)  =  V95  ttyrr1 


i/e  can  now  find  all  first-order  coefficients  u,  a  and  a  as 
functions  of  0« 

To  obtain  second-order  coefficients  we  first  solve  th8 
three  simultaneous  equations  for  U2?  a 2  and  o2 


-2u0a0 

u2  -  d  2  -  - 2_2 — _ 

*  d  ti«  -  an  -  1 


(3.5) 


(u2  +  a2  +  1)  +  (1  +  X)(u2  +  d2  ^  =  2(ao  +  ^u0) 

2 

a2  *  \d2  + 


where 


^  =  d£ 


(3.6) 

(3.7) 


The  results  are 


a0  -  ijiof  +  4\u0j 


“z  =  Tfi  VETf 
a2  -  5trVrr{a° '  3u° '  +  Xuo] +  itl°i 

°2  =  friVrrK  "  3u°  ‘  ili0i +  Xtt°) 

Solving  equation  (2,7)  gives  us 


[  u2e^  +  0  2e^®^  ** 


Cu2e(«)  +  °2e  (° )  ^ 

-fg(0)d9  rQ  j9  g(0)d0 
ea  h(e)e0 


d9 


g(0) 


=  fcos  e  4  (uoe  +  aoe)aln  el> 
[sin  0  -  (uo0  +  aoe)cos  e 


where 
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h(e)  = 


-  (u3e+ale)  f^ule+<3le^sln  9  + 


oeaoe 


sin  ©  =  (uoe+aoe)cos  © 


Solving  equation  (2.6)  gives  us 


U2e  “d2e  =  ^ 


(u0«-a0J  -  (tt^-a1a)eos  ©  £(uie“dle)sin  8 


ule**ale  2e  20 


‘le  le 


+  2(uie-dle)cos  8]  "  4uoeaoecos  8}  (3o9) 


With  the  equation 


p.d 


«2e  +«  2e  =  -  (£$(u2a  -  °2e^  +  l?\(u2e  *  a2e>  + 

(3.10) 

we  then  have  three  simultaneous  equations  to  solve  u2e?  2q 
and  a2e  explicitly  as  functions  of  ©. 

RftKlPB-Q 

The  equation  of  state  in  region  c  is 


p  "  l 

pd  [ 


where  n,  B  and  pQ  are  constants  determined  by  Richardson, 
Arons  and  ^alverson  (1947) 

n  =  7 

B  *  3.311  k  10  9  dyr./cra2 
P0  =  1.007  g/cm3 
pd  =  24.746  x  10 10  dyn/cm2. 


Equation  (3.11)  provides  us  with  the  corresponding  values 


of  Lc,  Xc  etc 
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The  computation  involved  in  solving  uc  and  ac  :is'  much 
simpler  than  those  in  region  e0  The  required  integration 
constants  and  boundary  conditions  are  supplied  by  Holt  (1956) 
and  coefficients  of  first  and  second-order  are  evaluated 
according  to  the  expressions  in  II. 

S£glPJL£» 

Coefficients  in  region  g  are  evaluated  in  a  similar  manner 
to  those  in  region  c  using  the  expressions  listed  in  II,  Since 
this  is  a  uniform  region  as  far  as  certain  properties  are  con¬ 
cerned,  only  first  and  second-order  coefficients  are  functions 
of  9  and  values  of  L  and  P  are  constants  equal  to  the 
values  of  XQ,  LQ  and  Pe  at  the  boundary,, 

All  values  of  uQ,  ae,  d0,  uc,  ug  and  ag  are  tabulated  as 
functions  of  9  in  V.  All  necessary  data  for  the  computation 
obtained  from  Holt  (1956)  are  given  in  the  appendix,, 

IV,  The  construction  of  characteristic  lines, 

The-mlnus  characteristic  line. 

A  fan  of  minus  characteristic  lines  goes  through  O’  in 
the  expansion  region,  among  them  Is  the  boundary  between  regions 
e  and  d.  The  boundary  Is  expressed  by  the  equation 

9  =  2.5997  +  Oo 61975  -  1.8084 +  ... 

Our  problem  now  Is  to  determine  the  fluid  properties  along  this 
line. 

We  know,  from  the  result  of  Coring,  that  in  the  detonation 
region,  the  dependent  variables  are  expressed  in  the  form 
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u  =  u0  +  u-^  +  UgC2  +  •  , 


Substitute  them  into  the  equations 

(u  +  a  -  z)  SLmjLaI  =  7  2U a  (4a  1) 

dz  z 

2 

where  z  =  1  ~  » 

V/e  can  determine  a  point  of  known  properties  by  choosing 
a  small  value  of  Z  ,  thus  of  @  and  £  .  From  this  point  on  we  can 
determine  all  the  properties  along  the  bounding  characteristic 
by  integrating  (4,1)  numerically. 

The  initial  point  is  chosen  at  (=  0,1  or  z  =  0,99,  u,  a 
and  d  are  evaluated  respectively, 

Equations  (4,1)  are  replaced  by  simple  differential  equa¬ 
tions  (4.2)  as  a  first  approximation 


— 2— +-d  2)  ~  (ul  +  d!)  =  _  1 

z2  "  zl  Z1  +  ax  -  zx 


(uo  -  d9)  -  (u^  -  o^) 


(4.2) 


2uia- 


z2  -  Z1  ux  -  ax  - 

z2  is  a  value  chosen  such  that  rapid  convergence  ensued. 

U2  and  d2  can  be  determined  from  (4.2)  and  a2  may  be  obtain¬ 
ed  from  the  relation 


0  =  i\  & 

k  pa 


A  better  approximation  of  u2  and  d  2  can  4®  achieved  by 
writing  the  difference  equation  in  the  form 
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(ug+ggMtil+dj)  ^  J  2uiai  +  2u2a2 

Z2“zl  " l ( u^/"” 2^z i^  z2( iig"^ a 2** z 2^ 


(4.3) 


(q2-g2)-(urd1)  _  J  2u1a1  g  e  f 

ZP“Z1  Vzl(ul“®l“zl^  z2(u2+a2’z2)  J 


2u2a2 


Values  of  u2,  a2  from  the  first  approximation  are  used? 
and  the  new  value  a2  may  be  found  by  using  the  new  value  of  Og® 
This  process  may  be  carried  out  a  few  more  times  to  obtain 
enough  accuracy,,  The  values  of  u,  a  and  c  at  different  z  or  C 
are  tabulated  in  Table  4,  V, 

The  plus  characteristic  line 

The  "plus"  characteristic  family  lie  across  the  flow  field 
from  the  boundary  between  e  and  d  to  the  main  shock,,  One  member 
of  the  family  lying  very  near  0f  is  constructed  in  the  follow¬ 
ing  manner o 

The  governing  equation 


sSL  =  £.5  ll-  ±  C ii . Jan, jb) 

d9  u  +  a  -  tan  6 


(4.4) 


is  replaced  by  a  simple  difference  equation.  A  small  value  of 
K  on  the  boundary  between  e  and  g  is  chosen  to  be  used,  as  an 
initial  point.  Having  found  all  the  coefficients  of  series 
expansion  in  regions  e,  g  and  c,  we  are  able  to  integrate  the 
above  equation  numerically  in  the  same  manner  as  the  minus 
characteristics.  This  process  gives  the  position  and  condition 
of  this  particular  "plus"  characteristic.  These  data  are  given 
in  Table  5,  V. 
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-Characteristic  lines 
Having  determined  one  "plus"  and  one  "minus"  characteristic 
line,  we  are  in  a  position  to  construct  the  whole  flow  field 
by  the  method  of  characteristics. 

Due  to  limit  of  time  the  whole  work  is  not  tried  at  the 
present  stage.  However,  a  plus  characteristic  line,  has  been 
computed  in-  region  e  by  IBM  CPC  as  an  illustration  of  the 
method. 

Using  a  modified  method  of  Berry,  Butler  and  Holt  (1955) 
we  introduce  new  independent  variables  x  and  y  by  the  relations 


x  =  t  -  1 
y  =  r  -  1 

The  characteristic  lines  are  defined  by 

4^  =  u  ±  a  (4,5) 

dx 

The  characteristic  equations  are  then 

du  +  dd  =  +  dx  (4,6) 

* *  y 

In  numerical  work  equations  (4,5)  and  (4,6)  are  replaced 
by  simple  difference  equations  and  an  iterative  method  may  be 
employed  as  in  the  construction  of  the  "minus"  and  "plus" 
characteristics. 

Figure  2  shows  a  characteristic  quadrangle  in  which  12  and 
13  are  given  segments  of  known  "minus"  and  "plus"  characteristics 
respectively  and  tho  point  4  is  to  be  constructed,  A  first 
approximation  is  made  by  calculating  the  position  of  point  4 
as  the  intersection  of  tangents  to  characteristics  at  points  2 
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and  3,  using  equation  (4,5)?  then  finding  values  of  u  and  d  at 
4  by  solving  equations  (4.6)  as  simple  difference  equations# 
Better  approximations  may  be  obtained  by  iteration  in  the  same 
manner  as  described  above,  The  whole  process  is  programmed  on 
an  IBM  CPC  computer  in  a  continuous  manner.  A  point  lying  out  .on 
the  minus  characteristic  was  picked  to  be  point  2,  all  the  needed 
known  properties  associated  with  this  point  are  stored  in  the 
machine,  likewise  are  the  data  of  a  chosen  point  3  on  the  known 
plus  characteristics.  The  two  difference  equations  (4. 5)  are 
solved  to  give  the  values  of  x,  y  and  then  u  and  d  in  (4.6), 

The  value  of  d  is  then  used  to  calculate  its  corresponding  value 
of  p  from  the  relation 


&  *  -JL 

dp  kpa 


(4.7) 


Equation  (4,7)  is  replaced  by  a  simple  difference  equation 
using  either  point  2  or  3  as  one  end  point.  The  value  of  p 
is  iterated  according  to  (4,7)  and  its  final  value  is  used  to 
calculate  a.  After  the  positions  and  properties  of  point  4  are 
found,  the  machine  automatically  stores  everything  connecting 
4  into  those  previously  containing  2  and  a  set  of  data  of  a  new 
point  3  is  fed  in  and  the  process  is  repeated,  A  line  contain¬ 
ing  roughly  25  points  may  be  constructed  in  a  matter  of  six 
hours.  It  is  thus  felt  that  CPC  is  not  the  most  suitable  equip¬ 
ment  for  doing  this  type  of  work,  both  time  and  costwise.  The 
value  of  the  "plus"  characteristic  is  given  in  Table  6,  V, 
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V*  Tables 


Table  le 

Table  2* 

Table  3* 

Table  4. 
Table  $e 

Table  6„ 


Coefficients  of  series  expansion  In  region  e,  u, 
a  and  d . 

Coefficients  of  series  expansion  In  region  g,  u 
and  a* 

Coefficients  of  series  expansion  In  region  c,  u 
and  a. 

Second  plus  characteristic  in  region  e« 
Arbitrary  plus  characteristic,  x,  y,  u,  a  and  d 
near  O’  throughout  the  field* 

Minus  characteristic,  x,  y,  u,  a  and  d  along 
boundary  between  e  and  d« 
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Table  1.  Coefficients  of  series  expansion  in  region  e 


u 


9 

u0e 

*le 

u2e 

m 

u0e 

ule 

u2e 

2.5997 

.1990 

-.3671 

-.023467 

2.81  , 

.2705 

-.2431 

,114960 

2.6o 

.1990 

-.3671 

-.028397 

2.82 

.2759 

-.2347 

.129732 

2,61 

.2012 

“.3635 

-.026864 

2.83 

.2817 

-.2272 

0 144403 

2.62 

.2032 

-.3599 

-.024931 

2.34 

.2875 

-.2188 

.161040 

2.63 

.2054 

-.3559 

-.022719 

2.85 

.2941 

-.2106 

.178349 

2a64 

.0076 

“.3518 

-.019769 

2.86 

.3005 

-.2028 

.195608 

2.65 

.2101 

-.3475 

-.016512 

2.87 

.3067 

-.1945 

.214380 

2.66 

.2127 

3428 

-.012531 

i 

2.38 

.3138 

-.1873 

.231502 

2.67 

.2153 

-.3382 

3 

iH 

CC 

O 

O 

0 

1 

2.8876 

.31915 

-.1819 

.250759 

2.68 

.2181 

-.3329 

-.003510 

2.69 

.2211 

-.3275 

+.001988 

2.70 

.2242 

-.3220 

+.007069 

2.71  j 

.2275 

-.3159 

+.014615 

2,72 

.2309 

-.3100 

+.020684 

2.73 

.2345 

-.3033 

+.028098 

2.74 

.2384 

-.2966 

.037045 

2.75 

.2424 

-.2903 

.045066 

2.76 

.2465 

-.2834 

.054294 

2.77 

.2508 

-.2763 

.064760 

2.78 

.2555 

-.2691 

.075793 

2.79 

.2602 

-.2611 

.087235 

2.80 

.2654 

-.2516  ! 

.100888 
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Table  I  (continued) 
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J 

,  | 

'i; 

■i 

I 

i 

t 

9 

a 

a 

a 

9 

a. 

i 

Oe 

le 

2e 

Oe 

lej 

2e 

■t 

\ 

2® *997 

.solo 

-.7893 

.600824 

2.81 

.61*2 

-.2488 

.36240 

t 

2»60 

.8010 

-.7893 

o*97736 

2.82 

.6096 

-.2323 

.37880 

1 

i 

2.6l 

.7890 

-.7481 

.**900 

2.83 

.6o4* 

-.2186 

.38830 

|  : 

2*62 

.7778 

-.7103 

.*114* 

2.84 

.*99* 

-.2033 

.40273 

,  i 
' ! 

2.63 

.7666 

-.'67*0 

.47830 

2.8* 

.*947 

-« 1894 

.41910 

1 

: 

2.64 

.7*60 

-.6414 

.44*00 

2.86 

.*900 

-.1767 

.43800 

■ 

2.6* 

.74*8 

-.6107 

.42480 

2.87 

.*8*0 

-.163* 

.4*8*0 

. 

' 

2.66 

.73** 

-.*703 

.4o46o 

2.88 

.*810 

-.1*22 

.47892 

1 

2.6  7 

.72** 

-.**12 

.38798 

2.3a76 

.*787 

-.1440 

.49168  j  ‘  1 

1 

2„68 

.71*9 

-.*226 

.37364 

: 

2.69 

.7065 

-.4954 

.360*0 

2.70 

.6070 

-.4709 

.349*0 

2.71 

.6886 

-„4i+?9 

.341*0 

! 

2.72 

.6800 

-.4234 

.314*0 

i  | 

2.73 

.6713 

~«4oo4 

.32800 

2.74 

.6632 

-.3790 

.324*0 

2.7* 

.65** 

-.360* 

.319*0 

2.76 

.6480 

-.3414 

.319*0 

*  - 

2.77 

.6410 

-.3229 

.32400 

2.78 

.6340 

-.30*1 

.331*0 

' 

2.79 

0627* 

-.2866 

.34200 

2.80 

.6210 

-.26*9 

.3*1*0 
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2  7 


\  Table 

r 

5 

i 


e 

d0a 

6  le 

2,5997 

1.612 

-.1294 

2.60 

1.6118 

-.1294 

'  2.61 

1.6096 

-.1324 

f  2.62 

1.6075 

-.1353 

|  2.63 

1.6052 

-.1386 

t  2.64 

i' 

1.6030 

-.1419 

r  2.65 

r 

I.6008 

-.1454 

f  2.66 

i 

1.5982 

-.1493 

•  2.67 

1. 5952 

-.1531 

[  2.68 

1-5925 

-.1574 

2,69 

1.5396 

-.1619 

■  '  2.70 

1.5865 

-.1664 

2.71 

1,5830 

-.1715 

|  2.72 

1.5300 

-.1764 

:  2.73 

1.5764 

-.1820 

2.74 

1,5727 

-.1876 

*  2.75 

1.5690 

-.1923 

I  2.76 

1.5650 

-.1985 

[  2.77 

1.5605 

-.2044 

l  2.78 

1,5560 

-,2io4 

\  2.79 

1,5510 

-.2171 

2.30 

1.5456 

-,2253 

I  (continued) 

6  Q 


4 


5 

5 

'6 


20  82 


15431  2 8 84  le5230  -,2526  -.4 

15744  2.85  1,5170  -02594  -.4 

6133  2.86  1.5104  -.2657  -.5i 

6556  2.87  1.5038  -.2725  -.5 

7057  2.88  1.4970  •7,2782  -.5! 


nlWn 
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Table  2,  Coefficients  of  series  expansion  in  region  g 


28 


u 


6 

°g 

Ulg 

2g 

@ 

uog 

U-. 

lg 

U2g 

2.8876 

.3191 

-.1881348 

,0321892 

3.24 

.3191 

-.140979 

.167320 

2.8880 

.3191 

-.18753696 

•%  0014837 

3.26 

.3191 

-.138320 

. 174408 

2.8885 

.3191 

-.18720721 

.0099878 

3.28  ; 

.3191 

-.135612 

. 181429 

2.8890 

.3191 

!  -.18695451 

.0146487' 

i  1 

3.30 

.31911 

-.132855 

.188381 

2.890 

.3191 

!  -.18654632 

1 

.0193006! 

3.32 

.31911 

-.130044 

.195259 

2.900 

.3191 

-.134051 

.0309208 

3.34 

.3191; 

1 

-.127178 

0 202060 

2.920 

.3191 

-.130703 

.0410614 

3.36 

.3191 

124252 

.208783 

2.9^0 

.3191 

| -.177898 

.0496793 

3.38 

.3191 

-.121263 

.215423 

2.960 

.3191 

-.175308 

.0579384 

3.4o 

.3191 

-.118209 

.221978 

2.980 

.3191 

| -.172829 

.0660172 

3.42 

.3191 

-.115083 

.228445 

3.000 

.3191! 

!-. 170410 

! 

.0739900 

3.44 

.3191 

-.111884 

.234819 

3.02 

.3191 

-.168025 

.0818851 

3.44051 

.3191 

-.110173 

.238129 

3.04 

.3191 

-,165653 

.0897163 

3.06 

.3191 

-.163233 

.0974877 

3.08 

.3191 

| 

-.160906 

.10865214 

3.10 

.3191 

-.158514 

.11613660 

3.12 

.3191 

-.156103 

.12358839 

3.14 

.3191 

-.153666 

.13100068 

3.16 

.3191 

-.151200 

.13836904 

3.18 

.3191 

-.148702 

.14568949 

3.20 

.3191 

-.146167 

.15295665 

3.22 

.3191 

-.143594 

.16016830 

56207/9 


Table  2  (continued) 


29 


£1 


9 

a 

“Og 

a. 

lg 

a2g 

9 

1 

aog 

a!g 

a2g 

2*3876 

.5787 

-.0959742 

'  -.-24267490 

3.24 

.5787 

-■.089799 

-.330294 

2*8380 

.5787 

-.09625375 

-.21407971 

3.26 

.  57Q  7 

-.088921 

-.333795 

2*8885’ 

.5787 

-.09639008 

-.22933902 

3.28 

.5787 

-.088003 

-.337153 

2,8890 

.578 7 

-.09647710 

-.23234037 

3.30 

.5787 

087045  j 

-.340371 

2.890 

.5787 

-.09663147 

-.23540690 

3.32 

.5787 

-o036046 

-.343445 

2*900 

.5787 

-.097259 

-.24372671 

3.34 

*5787 

-o035005 

-.346374 

2.920 

.5787| 

-.097643 

-.25143223 

3.36 

.5787 

-.083921 

-.349157 

2.94 

.5787 

-,007701 

-.25783400 

3.38 

.5787 

-.082793 

-.351792 

2,96 

.5787 

-.097602 

-.26330211 

3.40 

.5787 

-.081620 

1  -.354281 

2.98 

.5787 

-.097398 

-,26949212 

3.42 

.5787 

-.080400 

-.356619 

3-00 

.5787 

-.097113 

-.27497359 

3.44 

.5787 

-.079132 

-.358808 

3o02 

.5787 

-.096763 

-.28027525 

3.4505 

.5787 

i  —0  078446 

-.359897 

3e04 

.5737 

-.096354 

-.2854  13  52 

3.06 

.5787 

-.095394 

-.29039747' 

3.08 

.5787 

-.095365 

-,29736800 

3.10 

.5787 

-.094831 

-,30194900 

3.12 

.5787 

-.09235 

-.30640000 

3.14 

.5787 

-.093596 

-.31071900 

3.16 

.5787 

-.092916 

l 

-.31490500 

3.18 

.5787 

-.092196 

-.31895800 

3.20 

.5797 

1 -.091437 

-.32287500 

3.22 

.5787 

-.090633 

-.32665300 
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Table  3*  Coefficients  of  series  expansion  in  region  c 


30 


u 

a 

e 

— 

u0c 

ulc 

u2c 

9 

a0e 

a_ 

1c 

a2c 

3.4505 

.3191 

-.108258 

-.002841 

3,4505 

.9971 

-.004631 

-.012636 

3.46 

.3191 

-.107620 

-.006328 

3.46 

«9°71 

-.004595 

-.012609 

3.48 

.3191 

-.106246 

-.013667 

3.48 

.9971 

-.004516 

-.012546 

3.50 

.3191 

-.104829 

-.021001 

3.50 

.9971 

-.004435 

-.012479 

3.52 

.3191 

-.103365 

-.028321 

3.52 

.9071 

-.004352 

-.012407 

3.54 

.3191 

-.101855 

-.035641 

3.54 

«9q71 

-.0042 67 

-.012329 

3.56 

.3191 

-.100297 

-.042940 

3.56 

.9971 

-.004180 

-.012247 

3.58 

.3191 

-.098688 

-.050223 

3.58 

.9971 

-.004090 

-.012159 

3.60 

.3191 

-.097026 

-.057485 

3.60 

.0071 

-.003998 

-.012066 

3.62 

.3191 

-.095310 

-.064724 

3.62 

.9Q71 

-.003903 

-.011968 

3.64 

.3191 

-.003537 

-.071936 

3.64 

.9971 

-.003806 

-.011865 

3.66 

.3191 

-.091703 

-.079120 

3.66 

.9971 

-.003705 

-.011758 

3.68 

.3191 

-.080805 

-,086272 

3.68 

.9971 

-.003602 

-.011645 

3.70 

.3191 

-.087841 

-.093389 

3.70 

.9971 

*000^95 

-.011527 

3.72 

.3191 

-.085804 

-.100468 

3.72 

.9971 

-.003385 

-.011405 

3.74 

.3191 

-.083692 

-.107506 

3.74 

.9971 

-,003271 

-,011277 

3.76 

.3191 

-.081495 

-.114500 

3.76 

.9971 

-.003153 

-.011145 

3.78 

.3191 

-.079210 

-.121448 

3.78 

.9971 

-.003030 

-.011008 

3.7881 

.3191 

-*078258 

-.124248 

3*7881: 

.9n71 

-.002979 

-.010951 
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Tabl8  4.  Second  plus 


y 


X 

.00491874 

.0049530 

.00409898 

.'00503931 

.00508140 

.00512419 

.,00516652 

.00520091 

.00525443 

.00529001 

.00534324 

.00533979 

.00543675 

.00548351 

.00553005 

.00557925 

.00562902 

.00567888 


^.00272657 

-.00269038 
-.00265519 
-.00261963 
-.00258272 
-.00254540 
-.00250866 
-.00247119 
-.00243  389 
-.00239470 
-.00235604 
-.00231731 
-,00227743 
-.00223782 
-,.00219769 
-.002156 88 
-.00211486 
-.00207277 


31 
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32 


Table  4  (continued) 


X 

y 

..... ... . ...  _  .— .  _ j 

u 

a 

d 

.00573133 

-.00202848 

.246718 

.59148 

1.525861 

.00578250 

-.00198  521 

.252399 

.593557  ] 

1.520165 

.00583420 

-.00194143 

.2581 7% 

.589306 

1.514371 

.00588693 

-.00189669 

.264273 

.585222 

1.508258 

.00593944 

-.00185203 

.270360 

.581503 

1.502155 

.00599500 

-.00180462 

.276905 

.577802 

1.495503 

.00605073 

-.00175688 

.283813 

.574329 

1.438666 

.00613273 

-.00168628 

,294272 

.569529 

1.478180 

• 

.00618400 

-.00164190 

.300613 

.566875 

1  1.471822 
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Table  5«  Position  of  an  arbitrary  plus 
oharaoteristic  near  0f 


33 


0 

K 

9 

5 

9 

K 

0 

K 

2,62 

.036675 

2,85 

.0387349 

3.24 

0,0466291 

3.68 

0.0698348 

2,63 

,0367317 

2.86 

.0388466 

3.26 

0,0472913 

3.70 

0.0716605 

2,64 

.0367915 

2,87  . 

.0389712 

3.28 

0.0479925 

3.72 

0.0736531 

2,6? 

,0368551 

2.88 

.0390937 

3,30 

0.0487359 

3.74 

0.0758384 

2,66 

.0369207 

2,3848 

.03916207 

3.32 

0.0495252 

3.76 

0.0782476 

2,67 

,0360.892 

2.89 

,03°2295 

3.34 

0.0503645 

3.78 

0.0809201 

2,68 

.0370607 

2,90 

.0393636 

3.36 

0.0512586 

2,69 

,0371350 

2.92 

.0396421 

3.38 

0.0522130 

2,70 

.0372124 

2.94 

.0399347 

3.4o 

0,0532337 

2,71 

.0372937 

2.96 

.0402420 

3.42 

0.0543281 

2,72 

.0373724 

2.98 

.0405647 

3,44 

0.0555045 

2,73 

.0374574 

3,00 

,0409037 

3,446 

0,0558721 

2,74 

,0375^65 

3.02 

.0412597 

3.46 

0.05660 

2,75 

,0376385 

3.04 

.0416336 

3.48 

0.0574515 

2,76 

.0377334 

3.06  . 

,0420264 

3.50 

0.0533548 

2.77 

.0378313 

3.08 

0,0424392 

3.52 

0.0593143 

2,78 

.0379320 

3.10 

0v04 2 87 3 2 

3.54 

0,0603354 

2,79 

.0380356 

3.12  . 

• 

0.0433296 

3.56 

0,o6l4238 

2,80 

.0381421 

3.14  ; 

0.0438098 

3.58 

0,0625864 

2,81 

,0382422 

3.16 

0.0443155 

3,60 

0.0638311 

2,82 

.0383535 

3.18 

0.0448482 

3.62 

O.0651667 

2.83 

.0384683 

3,20  j 

0.04  54100 

3.64 

o. 0666039 

2,84 

.0385350 

3.22 

0.0460028 

3.66 

0. 0681551 
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Table  6»  Minus  characteristic 


X 

y 

u 

a 

d 

,00116620 

-.00067019 

01899 

.7524 

1.6025 

,00402125 

-.00267994 

,1800 

.7162 

1.5928 

,00960330 

-.00525267 

.1725  | 

.6888 

1.5832 

.0172790 

-.00925181 

.1639 

,6610 

1.5712 

0 0270417 

-.0143749 

.1542 

.6400 

1.5590 

,104652 

-.063729 

;  .1045 

,5810 

1.4970 

.193669 

-.158140 

00527 

.5534 

1.4300 
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AREP.qflAX 

(From  Holt  (1956)  ) 


Table 

1*  Initial  values  of  dependent  variables 

air 

water 

o 

a 

II 

Pflg  =  W65 

2o763xlO~3 

6.132xlO_1 

II 

o 

o 

3 

u0g  "  a0ea) 

9.219X10"1 

3  o  191X10"*1 

II 

60 

o 

<d 

a0e  ^  ^ 

Q<,989xlO"'2 

5.787X10"1 

a0c 

2.957X10"1 

9. 971x10**  ^ 

Table  2.  Constants  of  integration 


air 

water 

air 

water 

Cx  -l.738xlO“3 

-5« 504x10" 1 

c4 

~3.342xlO~3 

-9.457X10”1 

C2  -3 o 797x10 "4 

-5o291x10“2 

c5 

9«428x10"3 

9.619x1c"1 

C3  -6O433X10”1 

0 

C6 

60663 

9 

Gx  -7.243X10"3 

-'5.604X10'”1 

g4 

1. 274x10 ~2 

-9.056xl0“2 

G2  -9.580x10“3 

9.264xl0"2 

G 

5.444xl0“2 

-1.271 
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